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Abstract 



The 5 = 1/2 antiferromagnetic Heisenberg model on multi-leg ladders is investigated. 
Criticality of the ground-state transition is explored by means of finite-size scaling. The 
ladders with an even number of legs and those with an odd number of legs are distinguished 
clearly. In the former, the energy gap opens up as AE ~ J±, where J± is the strength of 
the antiferromagnetic inter-chain coupling. In the latter, the critical phase with the central 
charge c = 1 extends over the whole region of J± > 0. 

> 

o 1 Introduction 

^ Understanding of the ground-state criticality of one-dimensional quantum system has been enriched 
0> greatly in the past decade. One of the important breakthroughs was finite-size scaling based on the 
conformal field theory |T[ ; see for a review. Nowadays many attempts on generalization for higher- 
g dimensional systems are in progress. One of them is the study of coupled chains. The Heisenberg 
i ladder models are particularly of interest. The two-dimensional Heisenberg model might be explored 
as the limiting case of the ladder models. In addition, actual substances which realize the ladder 
O models have been developed experimentally: Sr 2n _ 2 Cu2 n 04 n _2, for example 0. 
The Hamiltonian of the 5 = 1/2 Heisenberg ladder models is given by 



TC = Hleg + Tuning, (1-1, 



where 



L—l n t 

'Hlcg = J £>x,y ' S x +l,y, (1-2) 

x=0 y=l 
L-lnt-1 

Tuning = J± £*x,y ' &x,y+l- (^-3) 

x=0 y=l 

In the present paper we treat the antiferromagnetic case, J > and J±_ > 0. Each component of the 
spins is defined by 

S" = i<r" (V = x,y,z), (1.4) 

where {c^} are the Pauli matrices. We define S ± = S x ±iS y . We impose periodic boundary conditions 
only in the x direction: 

S L , y = S 0t y for y = l,2,...,m. (1.5) 

Hereafter the system size L is even. 

In the present paper, we focus on the ground-state criticality of the ladder models. It was conjec- 
tured Q that the ladder models with even ni and those with odd n\ behave quite differently from each 
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other; for J± — 1, the former is massive, while the latter is massless. Here we discuss the criticality 
for arbitrary positive Jj_ and for general rij. 

It is well known that the antiferromagnetic Heisenberg chain {J±/J = 0) is massless |J. Many 
questions arise when we introduce the inter-chain coupling 7i rung : It has been quite controversial 
whether the system with ni = 2 becomes massive for any positive J±, or there exists a critical point 
at a finite value of J±; It has not been clarified whether the massless phase for odd n\ extends over 
the whole region of J± > 0, or a massive phase exists in the region < J± < 1; The universality of 
the massless phase has not been discussed either. 

The purpose of the present paper is to answer the above questions on the basis of finite-size 
scaling theory. We introduce our scaling ansatz in section |2.1| . The ansatz yields useful conclusions 
immediately. In particular, we conclude that the whole region of J± > is controlled by the stable 
fixed point at J±/J = oo. Thus we distinguish the systems with even ni and those with odd rti clearly: 
the former is massive for any J±, while the latter is critical with the central charge c = 1 for any J±. 
We present some results for even n\ and for odd n\ in sections 2J2 and |2.3| , respectively. We particularly 
show estimation of critical amplitudes. We give numerical confirmation of the scaling ansatz in section 
2~T|. In addition to that, we confirm the ansatz by means of perturbation theory in section BL 



2 Finite-size scaling of the energy gap 

In this section, we introduce the finite-size scaling form of the energy gap. We numerically confirm 
the scaling ansatz. We present our conclusions drawn from the ansatz. 

2.1 General arguments 

The central scaling ansatz of the present paper is written in the form 

AE(L) ~ -A(x) = J L ^- with x = —, (2.1) 
L x J 

where AE is the energy gap between the ground state and the first excited state, and A denotes the 
relevant scaling function. We show, in figure ??, numerical data obtained by the Lanczos method [[K] . 
The data are scaled well over a finite region of x (namely < x < 5) for n\ = 2, 4, and over the whole 
region of x for n\ = 3. We confirm the scaling form ( |2.1| ) by means of perturbation theory in section [3[ 
The ansatz ( |2.1|) implies that the inter-chain coupling is a relevant operator. The coupling 7i rung 
drives the system away from the point J±/J = 0, and makes the system renormalized to the limit 
J±/J — > oo as L — > oo; see figure ??. Thus the fixed point at J±/J = oo controls the whole phase of 
the region J± > 0. Apart from correction to scaling, the thermodynamic limit L — >• oo is equivalent 
to the limit J±/ J — » oo. 

We can thereby understand the difference between the systems with even ni and those with odd 
n\. For even m, the spins on each rung form an S = singlet in the limit J±/ J — > oo. The system at 



the fixed point J±/ J = oo is the S — chain with the energy gap of the order of J± For odd n ; , 
on the other hand, the spins on each rung form an S = 1/2 doublet in the limit J±/J — > oo. In the 
first-order perturbation of J, we have the same energy spectrum as of the S = 1/2 Heisenberg chain. 
(See section [O] for details.) Hence the system at the fixed point J±/ J = oo is the S = 1/2 Heisenberg 
chain, which is massless, or critical. Thus the whole phase of J± > for even (odd) n\ is controlled by 
the massive S = (massless S = 1/2) Heisenberg chain. 

In order to argue asymptotic behavior of the scaling function A(x), we can exploit the two limits 
J±/J —>■ oo and L — > oo. In the thermodynamic limit L — > oo, we must have non-divergent value of 



(b) 



Figure 1: scale 
Scaling plot of numerical data for the J±_ de- 
pendence of the energy gap. The abscissa is 
scaled LJ±/J, while the ordinate is 

scaled as LAE/J: (a) m eg = 2, 4 < L < 12; 
(b) n leg = 3, 4 < L < 8; (c) ni eg = 4, 
4 < L < 6. 



(c) 

the limit x — > oo. We then postulate that the asymptotic form differs for even n; and for odd n\ as 
follows: 

a / \ { Ax for even n t 

Alx) ~ < , . , . as x — > oo. (2.2) 

I o tor odd n.; 

Here A and 6 are appropriate constants. We can derive this from the behavior of the energy gap 
in the limit J±/J — > oo. For even ni, a first excited state is created by exciting one of the S — 
singlets on the rungs. The energy gap to this state is proportional to J±: AE ~ Jj_A(x)/x oc J±. 
Therefore we should have A(x) ~ x. For odd ni, the energy gap for the first excited state is that of the 
5 = 1/2 Heisenberg chain, and hence is proportional to J/L: AE ~ (J/L)A(x) oc J/L. This yields 
A(x) ~ constant. The above postulate (|2.2|) is actually observed in the numerical data in figure ??. 



2.2 Ladders with an even number of legs 

According to the above ansatz Q2.2D , we have A(x)/x — ► A as x — > oo for even Hence in the 
thermodynamic limit L — > oo, the scaling form ( |2.1| ) is reduced to 

~ AJ ± ^ with i/ = 1. (2.3) 

The coefficient A defined in ( |2.2| ) turns out to be the critical amplitude. Thus the critical system at 
J±/J = (namely, a set of the independent Heisenberg chains) becomes massive immediately when 
we turn on the inter-chain coupling. The critical-point estimate J±/J = and the exponent estimate 



Figure 2: renom 
The conjectured renormalization-flow diagram 
for the Heisenberg ladders. The fixed point 
J±/J = is unstable with respect to the 
inter-chain coupling. The stable fixed point 
at J±/J = oo is described by the massive 
S = chain for even n leg , and by the mass- 
less S — 1/2 chain for odd n\ eg . We thereby 
conclude that the whole phase of Jj_ > is 
massive for even rii eg and is critical for odd 
ni e g. 



Let us estimate the critical amplitude A for rii — 2. We utilize the following general finite-size 



scaling form for one- dimensional quantum ground-state phase transitions [|12 

-Cx 



t \p(j\ \ Ax + D + e~ in the disordered phase, 

' ~ r> at the critical point, (2.4) 

in the ordered phase. 



J 




Here x = L\e\ v denotes the relevant scaling variable with e being the distance from the critical point; 
x = LJ±/J in the present case. The coefficient A gives the critical amplitude in ( |2.3| ). The coefficient 
C, on the other hand, gives the critical amplitude of the correlation length: ~ C(J±/ J) u . We 
introduced in [O] the amplitude relation 

A/C = v s , (2.5) 

where v s is the sound velocity. The coefficient B is predicted from the conformal field theory as 
B = v s 7tt], where 77 is Fisher's correlation exponent | 13[| . For the Hamiltonian ( |1.1|) with J±/J = 0, 
the Bethe-ansatz solution |14| gives 



hence we have 



7T 

77 = 1 and v s = -; (2.6) 
^ = ^ = 1.570796 ■• • and B = = 4.934802 ■•■ . (2.7) 



The coefficients D + and _D_ may have weak x dependence [|i2 



The asymptotic behavior in ( ^.2| ) for even n\ is consistent with the scaling form ( |2.4| ) in the disor- 



dered phase. We can employ our previous analysis in fl2"f , where we showed for the one-dimensional 



quantum Potts model that a good estimator of the amplitude A is given by 

AT - mm 1^11) . (2 . 8) 

(If the data followed the scaling form Q2.4p completely, the quantity LAE/(Jx) would converge to A 
exponentially as x — ► 00. However, there appears the minimum in practice, probably because some 
correction to scaling fll2| .) The estimator (|2.8| ) may converge to A exponentially |l^] in the form 



AT ~ A + ci exp(-c 2 L). (2.9) 

We calculated A™ in for the data in figure ?? (a), and fitted the results to the form ( [2.9|) ; see figure ??. 
We thereby have the estimate A = 0.47(1) and hence 1/C — v s /A = 3.34(7) for n\ = 2. (The error of 



(a) (b) 

Figure 3: fit 

Estimation of the critical amplitude A for n\ eg = 2 as was done in ||12|| . (a) The quantity LAE/(Jx) 
versus x for L = 12, for example. The estimator (|2.8| ) is obtained as A™ n = 0.505 • • •, which is 
indicated by the dotted lines, (b) The values of the estimator Q2.8p are fitted to the form (|2.9p. The 
solid line shows the best fit Af in ~ 0.468 + 0.505 exp(-0.218L). 

the estimate A was evaluated through the least-squares fitting to the form (|2.9|).) Namely, we have 

AE~ 0.47 Ji and £-3.34^- (2.10) 

as L — > oo and near J±/J ~ 0. These explain the following estimates in |3J quite well: AE = 0.504J 
and £ = 3.19(1) for rii = 2 and J±/J = 1. Although it is unnecessary that the critical behavior (|2.10| ) 



near J±/ J ~ is observed even for J±/J= 1, we see in figure 2 of that correction to ( |2.10| ) may 
be quite small in the region J±/J< 1. 

Incidentally, the energy gap is AE = J_l for rii = 2 in the limit J±/J — > oo, where the system is 
reduced to a set of independent dimers. A crossover from the behavior (|2.10| ) to the behavior AE = J± 
may occur in the region J±/J> 1. 

If we naively apply the same analysis as above to the data for n\ = 4 with L = 2, 4 and 6, we 
obtain the tentative estimates 

A ~ 0.27 and C' 1 = v s /A ~ 5.8. (2.11) 



Though there may be some errors in these estimates, the values ( 2.10|) and ( p. 11] ) are fairly consistent 
with the scaling hypothesis £ oc Q . 

2.3 Ladders with an odd number of legs 

For odd ni, we can see in figure ?? that the scaling region is quite wide. We may naturally assume 
that the scaling form ( |2.1| ) with (|2.2| ) is valid for any value of L and J±/ J. Thus we have the following 
remarkable conclusion: for any positive J±, the energy gap in the thermodynamic limit L — > oo behaves 
as 

AE b , 

— s !• (2 - 12) 

where the amplitude b is the coefficient defined in (|2.2|). Note that this amplitude is independent of 



J_i_ is identical to that for J±/J — > oo, or the spectrum of the 5 = 1/2 antiferromagnetic Heisenberg 
chain. We thus conclude that the whole phase of Jj_ > for odd rii is a critical phase with the central 
charge c = 1. Note that the central charge is c = rii for J±/J = 0, because we have rii systems of c = 1 
at this point. 

We can obtain the value of the amplitude b by considering the limit J±/ J — ► oo. For this purpose, 
we first describe how to obtain explicitly the effective Hamiltonian of the 5 = 1/2 Heisenberg chain in 
the limit J±/J^ oo. 

In the very limit of J±/J — oo, or J/J± = 0, the Hamiltonian is reduced to H Iung ; the rungs are 
independent of each other. The ground state of the whole system is the direct product of the ground 
state of each rung. On each rung, an odd number of spins are coupled with the antiferromagnetic 
interaction J±. Hence two ground states of each rung are degenerate with 5=1/2 and S z = ±1/2. 
We express the states on the xth rung as \+ x ) and \— x ). The ground states of the whole system have 
the 2 L -fold degeneracy. This degeneracy is lift up in the first-order perturbation of J/J±, or H\ eg . We 
calculate the first-order perturbation for degenerate states by writing down the secular equation: 



H — XI 



0. (2.13) 



Here the matrix J denotes the identity operator. The first-order energy is denoted by A. The matrix 
H is a 2 L x 2 L matrix which represents the operator H\ eg in the subspace of the degenerate ground 
states. For example, one of the diagonal elements of H is given by 

( — 1 — 2 +3 — 4 ' ' '| K-lcg \~1 ~2 +3 — 4 ' ' ') 

L-l rii 

— J /, /j ( — 1 — 2 +3 ~4 " " "| S Xj y ■ S x+ i t y | — i —2 +3 — 4 ■ • •) 

x=0 y=l 

= JY, (<"il hi) HI S ly I"*) + HI S lv (+sl S lv l+8> + •••)' ( 2 ' 14 ) 

y=l 

and one of the nonzero off-diagonal elements is given by 

J n ' 

(-i -2 +3 -4 ■ ■ -I n lcg h + 2 -3 -4 ■ • •) = - E (-2I s Vy 1+2) (+ 3 | sl y 1-3) . (2.15) 

/ y=l 

It is apparent that the secular equation ( 2.13|) is equivalent to the eigenvalue equation for the 
5=1/2 Heisenberg chain: 

L-l 

7~(-cS = J c r X/ ■ S x+ \. (2-16) 

x=0 

We can see in fl2.14|) that the effective coupling J e s is given by 

nt / n2 

4 = je((+i 5 ;i+)) • (2-i7) 

y=l 

We show numerical results for J e g in figure ??. The value of J e s is generally close to unity but increases 
monotonically. In particular, J c g = J for n\ = 3 [11]. For large ni, the matrix element in (|2.17 ) may 
depend on n\ as M 

( + \S z y \+)~nr? /u = nr 1/2 . (2.18) 

Hence we may have J c ^ ~ 0(ni°). 

As we mentioned in ( |2.4| ) and (277), it is known [TJ] for the antiferromagnetic Heisenberg chain 
with the coupling J ef j that the energy gap behaves as 



v.nri , n 2 J, 



off' 



Figure 4: Jeff 
The effective coupling J cS for odd ny eg in the 
limit Jj_/J — > oo. The values were calculated 
numerically on the basis of ( p. 17 ). 



Figure 5: spectrum 
A schematic view of the low-lying energy spec- 
trum of the S = 1/2 antiferromagnetic Heisen- 
berg chain. 



Comparing this with ( [2.12 ), we arrive at the conclusion 



b = lim A(a;) = v s 7iri— = — (2.20) 

3 Perturbational derivation of the scaling ansatz 

In this section, we describe a perturbational calculation which yields the scaling form fl2.1|) . This was 
briefly reported in ||. 

3.1 Zeroth order of J±/ J 

Let us first describe the energy spectrum for J±/J = 0. At this point, we have n\ number of the 
S = 1/2 Heisenberg chains independent of each other. The spectrum of each chain is given as 
figure ??. Let |sj,) denote the singlet ground state of the yth leg, and let \M y ) = {\l y ) , \0 y ) , |— l y )} 
denote the triplet of the first excited states. The energy of the singlet ground state is written in the 
form 

V IXC ( 7T^ \ 

E s ~ e Q LJ ^-j-J = I e L - Y^-jr J J as L -> oo, (3.1) 

where e = 1/4 — In 2 is the exact ground-state energy density 0, v s = n/2 is the sound velocity as 
before, and c = 1 denotes the central charge. The energy gap to the triplet states is given by 



where rj = 1 is Fisher's correlation exponent as before. 

The ground state of the whole system is given by the direct product of the singlet states as 



W,(o) 

S^gs 



\ s y) 



y=l 



with the energy 



The first excited states have the 3nrfold degeneracy: 

/ 



J. 



Vff) = \M yi ) 



\ y^vi 



\ s y) 



(M = 0, ±1 and y 1 = 1, 2, • • ■ , m). 



The energy gap to these excited states is the same as in 



(3.3) 



(3.4) 



(3.5) 



3.2 First-order perturbation of J±/ J 

Now we calculate the first-order perturbation with respect to Tt rung . The first-order energy of the 
ground state is given by 



gs \ V-gs 

because we have 



n 



rung 



L-Xni-1 

^ff) = J ^Y. £ ( s v\ &.v \ s v) ■ ( s y+±\ &*,v+i = °> 

x=0 y=l 



(s|^|s) = (s|5 x ± |s) = 0. 



(3.6) 



(3.7) 



In order to obtain the first-order energy of the excited state, we have to construct the secular 
equation, because the zeroth-order first excited states are degenerate. We concentrate on the sector 
J2 y M y = 1. The ground state of this sector gives the perturbed first excited state, because the 
perturbed energy spectrum would be similar to that in figure ?? owing to the SU(2) symmetry. The 
following rii states are of this sector: 



\(Vi)) = IO® 



(yi = 1,2, •■•,m). 



y=l 



We immediately have 



((yi) | Wrung I (1/2)) = unless \yx -y 2 \ = l. 
The nonzero matrix elements are 



(3.8) 



(3.9) 



J 1 



L-l 



((y + 1)| Wrung \(y)) = ( S V I S x,y\ l y) <Vn I S tv+l\ S V+l) = J ± a 

Z x=0 

( 3 / = l,2,...,n I -l) 

and their conjugates. Here the coefficient a is defined by 

L 



a = — 
9 



(s|^|l> 



(3.10) 



(3.11) 



(The x dependence of the matrix element (s| S~ |1) appears only in the phase factor in the form e lkx , 
and hence the absolute value is independent of x.) The secular equation for the first-order energy of 
the excited state is thereby written in the form 



-Eg J±a 
J ± a -£« J ± a 







J±a 







J±a 

J±a -Eg 



0. 



(3.12) 



This is equivalent to the one- dimensional tight-binding model under the free boundary condition; the 
triplet state hops from a leg to a neighboring leg with the hopping amplitude J±a. This model is 
exactly solvable [15|. The zeroth-order wave function and the first-order energy are obtained in the 
forms 



(0) 



2 ^ . 

L sin 

y=l 

—2J±a cos 



urn 
ni + V 

TflTT 



y)\(y)), 



ni + l 



(3.13) 
(3.14) 



for m = 1, 2, • • • , m. The state with m = 1 gives the lowest energy. 
We thereby obtain the energy gap up to the first order of J±/ J as 



AE 



J 
L 



7T 

T 



2xa cos ■ 



7T 



ni 



(3.15) 



where x = L J±/ J as before. 

We are in position to estimate the size dependence of the coefficient a. In the following we show 
analytically and numerically that a = 0(L°). 

We first show that the coefficient a is proportional to the transverse susceptibility of the Heisenberg 
chain. Let us consider the Hamiltonian with a transverse field at the origin: 



L-l 

AeH = J X; Si ■ S i+ i 

i=0 



After the standard calculation, we have 





Xoo 



7=0 

dr(^(r)> 



E 



(3.16) 



! 



(3.17) 
(3.18) 



where the angular brackets (• • •) denote the expectation value with respect to the ground state of 
( p. 16 ), AE n is the excitation energy of the excited state \^ n ), and 



S x (r) = e 



D e 



(3.19) 



The most dominant of the excited states \ip n ) is \S — 1, S z = ±1) with the energy gap ( p.2|) , or AE 
L^ 1 . Hence the expression ( |3.17 ) is approximately rewritten as 



Act 



Figure 6: coeff 
The size dependence of the quantity a appears 
to be of the form a ~ (lnL) 1 ^. The solid 
line connects the last two data points with the 
slope uj ~ 0.68. 



On the other hand, a field-theoretic description []T6| yields 

1 

'o ~ -■ 



(3.21) 



Thus the expression ( |3.18p is a dimensionless quantity: Xoo — 0(L°). We thus arrive at the conclusion 
a = 0(L°). There may be logarithmic corrections to this size dependence. 

We also calculated the coefficient a in ( p. 11 ) numerically for L < 24; see figure ??. The logarithmic 
plot reveals that the coefficient actually behaves as 



a ~ (hiL) 



0.68 



0(L C 



(3.22) 



3.3 Order estimation of higher-order perturbations 

We next derive the approximate second-order energy of the ground state. The second-order perturba- 
tion is given by the formula 

£{2) _ _ \yPn\ Tuning iVwj ^ 23) 

M^m) AEn 
The most dominant excited states are the 3(ni — 1) states 

f } ( 

\^ n ) = { |1 W -W> }®\ (g) Is,)! (3.24) 







( nl \ 












l y=i / 



with ?/! = 1, 2, • • • , rii — 1- The energy gap Ai? n to these states is twice as large as in ( |3.2| ): 

A£ n ~ (3.25) 

The numerator of ( p. 23 ) gives the factor a 2 for each of the states Q3.24 ). Hence the second-order energy 
of the ground state is 

Eg) ^j^K^J_ (3 26) 

Although it is complicated to calculate the second-order energy of the first excited state ( |3.13| ), we 
can see that the size dependence is the same as in fl3.26|) : 



/ LJ±\ 2 J 



We thereby have 

AE ~ — 
L 



, ,2 

— + a\ax + a 2 {ax) 



(3.28) 



where x = LJ±/J, and «i and «2 are constants. Since a is of the order of L°, the expression fl3.28|) is 
consistent with the scaling form (|2.1|) . 

It is generally difficult to calculate higher-order perturbation explicitly. However, it is possible to 
estimate the size dependence roughly |I7| . The A;th-order energy is approximately given by 

E(k) ~ H (V'gal ^rung (^2 | ^rung " ' " (0fc-l| ^rung |^gs) 

x[(AE 1 )(AE 2 )---(AE k _ 1 )}- 1 . (3.29) 

We may estimate the dimensionality of the operator TC run g at LJ± x [S] 2 ~ J±, because the magnetic 
operator S of the antiferromagnetic Heisenberg chain has the dimensionality [S] = L~^l v = L^ 1 ^ 2 
(This rough estimate is consistent with the estimation a = 0(L°).) On the other hand, each energy 
denominator AE is of the order of J/L, because the Heisenberg chain is critical. We thereby have the 
rough estimate 

E^ ~ J ± k (L = ) k . (3.30) 



,LJ L \ J 

Hence the energy gap may be given by the form 



AE^ij^^ (3.31) 



L k 



with appropriate coefficients This is consistent with the scaling ansatz ( |2.1| ). 

We have not confirmed that the above perturbational expansion is convergent. However, the 
numerical results in figure ?? do suggest that the series converges at least over a finite region of 
x = LJ±/ J. 



4 Summary 

In the present paper, we introduced the finite-size scaling form of the energy gap of the antiferromag- 
netic Heisenberg ladder models: 

AE{L)/J ~ L^AiLJ^/J). (4.1) 

We confirmed this scaling form numerically as well as by means of perturbation theory. On the basis 
of the scaling theory, we discussed the criticality of the ladder models in a unified way. The difference 
between the ladders with even ni and with odd rii was attributed to the different asymptotic behavior 
of the scaling function in the limit L J±/ J — > oo. 

For even ni, the energy gap develops in the form AE ~ J±. The whole region of J± > is a 
disordered phase with a unique ground state. (This ground state is reduced to a set of independent 
singlets in the limit Jj_/J — > oo.) We estimated for rii = 2 the critical amplitude of the energy gap 
and the correlation length around the critical point J±/J = 0. 

For odd ni, on the other hand, the whole region of J± > is the critical line. The energy 
gap vanishes in the form AE ~ L^ 1 for any J±. The critical line is controlled by the S = 1/2 
antiferromagnetic Heisenberg chain which is obtained in the limit Jj_/J — > oo. The central charge of 
the phase is hence unity. 
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